NMR evidence for a "generalized spin-Peierls transition" in the high magnetic field 

phase of the spin-ladder Cu 2 (C5N 2 Hi 2 ) 2 Cl4 
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The magnetic field-induced 3D ordered phase of the two-leg spin-ladder Cu2(CsN2Hi2)2Cl4 has 
been probed through measurements of 1 H NMR spectra and 1/Ti in the temperature range 70 mK 
- 1.2 K. The second order transition line T C (H) has been determined between H c \ = 7.52 T and 
H C 2 = 13 T and varies as (H — H c \) 2 ^ 3 close to H c i- From the observation of anomalous shifts and 
a crossover in 1/Ti above T c , the mechanism of the 3D transition is argued to be magnetoelastic, 
involving a displacement of the protons along the longitudinal exchange (Jm) path. 
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Two-leg 5=1/2 ladders are ID objects formed by 
two antiferromagnetically (AF) coupled Heisenberg spin 
chains. In zero external magnetic field, their ground state 
is a collective singlet state (5 = 0), separated by a gap A 
from the first excited states which are triplets (5 = 1) [Q. 
As a consequence, the spin-spin correlations remain of 
short range even when T — > 0, in spite of the strong inter- 
actions. There is currently considerable interest in these 
systems, often named spin-liquids, since the short-range 
singlet correlations of the ground state are believed to 
lead to superconducting correlations when mobile charges 
are added M. 

The fascinating properties of spin-liquids can also be 
revealed through the effect of a magnetic field H. This 
can be described in four steps: (1) For H ^ 0, the gap 
is reduced as A(H) = A — gfiBH. (2) At the so-called 
quantum critical point H — H c \ — A/g^sH, the spin- 
gap vanishes. At T = 0, this defines a (quantum) phase 
transition between gapped singlet and gapless magnetic 
phases. (3) For H > H c i, the gapless spin system still 
exhibits ID behavior at finite temperature, but as T is re- 
duced the magnetic correlation length and the spin-spin 
correlation functions now diverge (Luttinger liquid be- 
haviour) . This behaviour can be observed up to a satura- 
tion field H C 2 where all spins are polarized by H. (4) For 
H c x < H < H C 2, the transverse coupling J t between 
ladders should drive the system towards a 3D magnetic 
ordering at low T. The nature of the 3D phase, in the 
vicinity of the two quantum critical points, is expected 
to be highly unconventional. 

Points (1-3) were previously observed in NMR studies 
of the spin-ladder Cu^CsN^H^^CU || in which 
the low values of the AF exchange coupling (between 
spins 1/2 on Cu 2+ ions) along the legs (J|| ~ 3 K) and 
along the rungs (J± ~ 13 K), lead to experimentally ac- 
cessible values of H cl (~ 7.5 T) and H c2 (~ 13.5 T) [§. 
As to point (4) , specific heat measurements M j| in the 



field range 7-12 T have indeed revealed a phase transi- 
tion towards a 3D ordered phase for T < 1 K. However, 
no microscopic experimental insight has been reported 
so far, although this phase currently generates a large 
interest @0-|l|. 

In this Letter, we present a 1 H NMR study of 
Cu 2 (C 5 N2Hi2) 2 Cl4 in the field range 7.5-14 T, includ- 
ing the T-dependence (in the range 70 mK-1.2 K) of the 
lineshape and of the nuclear spin-lattice relaxation rate 
1/Ti. From the splitting of NMR lines, we define the 
transition line T C (H) below which 3D ordering occurs. 
In addition, we observe through 1/Ti a drastic change 
in the low-energy spin excitations below ~1.3 K, which 
is above T c . This behavior is correlated with anomalous 
shift of some 1 H lines, which we attribute to the displace- 
ment of protons involved in the exchange path along the 
legs of the ladder. This is argued to demonstrate the 
magneto-elastic nature of the transition, which is in some 
way analogous to the incommensurate magnetic phase of 
spin-Peierls systems. 

Experiments have been performed on a single crystal 
placed inside the mixing chamber of a 3 He- 4 He dilution 
refrigerator, the 6-axis of the crystal being parallel to H . 
In this orientation, the number of inequivalent 1 H sites 
in the crystal is reduced to 24. All sites experience dif- 
ferent hyperfine fields through their dipolar coupling to 
the electronic spins localized at the Cu sites fli|| . For 
large polarization of the electronic moments, these cou- 
plings lead to 1 H spectra extending over several MHz, 
which have been recorded at fixed field by sweeping the 
frequency. In Fig. [I] is shown the low frequency part 
of such spectra recorded at 70 mK at various values of 
H . One clearly observes a splitting of all lines starting 
at H = 7.55 T and increasing for higher H values. This 
is the signature of an ordered magnetic phase. Follow- 
ing the evolution of the spectrum with H for different T 
values allows an accurate determination of the transition 
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FIG. 1. Upper panel: low frequency side of 1 H spectra 
recorded at fixed values of H at T = 70 mK. The onset of 
the transition at H — 7.55 T is revealed by the splitting of 
the lines. Lower panel: low frequency side of 1 H spectrum at 
various T for H— 7.65 T; the transition occurs between 350 
and 245 mK. 

line T C (H). As shown in Fig. |^, T c rapidly increases as a 
function of H — H c \ and then saturates around 900 mK 
for H > 9 T. In this range, the transition was determined 
at fixed value of H and decreasing T. Again, we observe 
a line splitting (Fig. |l|), which quickly increases below T c 
and then saturates at lower T, as expected for an order 
parameter. The resulting experimental phase diagram is 
shown in Fig. |^. 

There are a few theoretical calculations of the line 
T C (H): Close to H c \. Giamarchi and Tsvelik predict 
a variation as (H — H c i) 2 ^ 3 , resulting from the con- 
densation of dilute hard core bosons JLOj. Wessel and 
Haas rather propose an (H — Hd) 1 / 2 variation As 
shown in the inset to Fig. 0, T c can be well-fitted to 
oc (H - i? c i) 2/3 , with H cl = 7.52 T. Note that the pre- 
diction of a first order transition with a jump of magne- 
tization for values of H close to H c \ [[ll) is not observed 
in our data even very close to H c \. In all cases, each 
line splits at T c keeping the same center of gravity, thus 
indicating that the magnetization is continuous with T 
or H , in agreement with earlier thermodynamic measure- 
ments Q. 

A careful examination of the whole 1 H spectrum 
(Fig. ||) above T c reveals that the T-dependence of the 
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FIG. 2. Phase diagram in the H — T plane as determined 
from NMR experiments. The shaded area corresponds to the 
T-range in which proton motion is observed (see text). The 
inset shows the variation of T C (H) in the vicinity of H c i, the 
solid line is a fit to oc (H — H cl ) 2/Z . 



shifts of the two lines at the left side of the spectrum 
does not scale with that of other 1 H lines. These two lines 
are assigned to the protons H(2) and H(4) ||, which are 
along the exchange path Jm corresponding to the atom 
sequence Cu-N-H- • -Cl-Cu (see inset to Fig. [l]). Since the 
shift of a proton 1 H(z) is given by Sh(i) = A(i)\cm m 
which A(i) is its hyperfme held and XCu the spin suscep- 
tibility per Cu atom, the absence of scaling can only be 
explained if A(2) and A(4) become T-dependent. This, in 
turn, can only occur if the distances H(2)-Cu and H(4)- 
Cu change jLj|. This is thus clear evidence that some 
kind of lattice instability occurs prior to the magnetic 
ordering. 

Since these protons are along the exchange path corre- 
sponding to J|| , any modification of the hydrogen bond 
should clearly change the magnetic excitation spectrum 
of the system, which can be probed by the nuclear spin- 
lattice relaxation rate (1/Ti). 

Such change is indeed observed in T\ data, measured 
between 1.2 K and 70 mK and reported in Fig. |]. There 
are three striking features in these data: i) the huge de- 
crease of 1/Ti up to 5 orders of magnitude at H— 8.0 T 
and 10.85 T. This decrease can be fitted by a power law 
(1/Ti oc T~ 5 ) Q. it) the increase of 1/Ti, attributed to 
the divergence of the spin correlation functions || , stops 
around 1.3 K for all field values between H c i and H C 2- 
Hi) 1 jT\ starts decreasing before the onset of the 3D tran- 
sition (as detected by the modification of the lineshape). 

We now discuss the possible nature of this 3D ordered 
ground state. From a theoretical point of view, the spin- 
ladder Hamiltonian can be transformed into an interact- 
ing spinless fermion model through the canonical Jordan- 
Wigner transformation [T?}] . In this representation, H 
acts as the chemical potential ^i, and for H = 7J c i, \i lies 
exactly at the bottom of the band. Increasing H further 
fills the band in. Since the value of the Fermi- wave vec- 
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FIG. 3. T-dependence of the shift K(T) of various proton 
lines (the full spectrum is shown in the inset), normalized 
at 6 K. K(T) scales to the uniform spin susceptibility for all 
lines, except for those corresponding to protons H(2) and H(4) 
which strongly deviate below T ~ 1.5 K. This can only be ex- 
plained if the position of these protons, hence their hyperfme 
field, varies in this T-range. 



tor fcp is set by the field, kp is incommensurate (IC) with 
the underlying lattice, except at half filling. Due to the 
divergence of the spin susceptibility at 2kp, the on-site 
magnetization of the ordered phase is also expected to 
be incommensurate. Between H c \ and -ff C 2, and at suf- 
ficiently low T, the low energy properties of the system 
are those of a Luttinger liquid | fl7|| . 

In the same field range, the spin-ladder Hamiltonian 
can also be approximately mapped onto that of an XX Z 
5=1/2 chain jl 12|. In this latter representation, an effec- 
tive spin 1/2 is introduced, whose eigenstates correspond 
to the singlet and the lowest state of the triplet on a rung, 
and the effective field H c g is equal to H — (H cl + H c2 )/2. 
In the following discussion, we shall use either the spin- 
less fermion or the XXZ language. 

It is well known that there are two possibilities to 
achieve 3D ordering at finite T for quantum spin chains: 
a transverse magnetic coupling J t , leading to some kind 
of AF order when J t £| ~ fc^T or a spin-Peierls (SP) tran- 



sition in presence of magneto-elastic coupling [|19[ . In the 
latter modulation of the lattice occurs, which is 

stabilized by the energy gain due to the opening of a gap 
in the magnetic excitation spectrum. The 3D character 
of the transition arise, in this case, from to the 3D nature 
of the elastic modes. 

The case of a transverse magnetic coupling for an as- 
sembly of ladders has been treated by Giamarchi and 
Tsvelik 10 1. In their model, the 3D ordering corre- 



sponds to a freezing of the XY degrees of freedom of 
the triplet states, and below T c , the local magnetization 
M z (R) is incommensurate along the direction of the lad- 
ders. A magneto-elastic scenario has been treated by 
Nagaosa and Murakami jflj, who considered a modula- 
tion of the exchange along the rungs Jj_ and by Calem- 
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FIG. 4. T-dependence of proton 1/Ti at different values 
of H . Left panel: data at H = 8.0 T and 10.85 T, measured 
on line II (see inset to Fig. 3 ). Data above 1.1 K are taken 
from M. Right panel: data at H — 7.65 T measured on line I. 



czuk et al. [|3| who found that a modulation of J» bet- 
ter explains specific heat data pj§ . As in the purely 
magnetic scenario, the local magnetization M Z (R) is IC 
along the ladder direction, and the 3D ordered phase 
is in some sense similar to the IC magnetic phase ob- 
served in spin-Peierls systems above their threshold field 
H c |19|. There is, however, a noticeable difference: in 
regular SP compounds, there is a commensurate (dimer- 
ized) phase, which is a collective singlet for < H < H c . 
In the spin-ladder system, the commensurability occurs 
for H = (H c i + H C 2)/2 (i.e. H c s — 0). Any extension 
of the commensurability around this H value would cor- 
respond to a plateau in the magnetization (which is not 
observed in Cu2(C5N2H 12 )2Ci4). Along the same line, 
the parts of the phase diagram close to H c i or H c2 in 
the ladder case correspond to a field range close to the 
saturation of the magnetization in the case of a regular 
SP system. 

We now compare our data to the predictions of these 
different models. All of them predict an IC modulation of 
M z (R) , giving rise to an infinite number of inequivalent 
sites. This should transform each NMR line into a double 
horned lineshape |p0|| . Because of the high density of 1 H 
lines in our spectra, we cannot distinguish whether each 
line transforms this way or simply splits. 

At variance with the lineshape, the T-dependence of 
1/Ti is expected to depend strongly on the nature of 
the ground state. A purely magnetic ground state [|t|Jic|] 
implies a divergence of 1/Ti at the transition, which is 
not observed experimentally. 

The increase of 1/Ti upon cooling, seen on Fig. 4 above 
~ 1.3 K, is indeed related to the Luttinger liquid behav- 
ior of the gapless ID system, as explained in ||] . This in- 
crease cannot be attributed to critical fluctuations linked 
to the transition, since it starts too far above T c Q and 
furthermore we now find that it even stops above T c . 
This is particularly obvious for the data at 7.65 T (right 



panel of Fig. 4). 

In contrast, in a spin-Peierls like transition, the low- 
energy spectral weight of AF fluctuations starts being 
suppressed even above the 3D ordering by the coupling 
to the elastic degrees of freedom. Hence, there is no di- 
vergence, but a rapid decrease of 1/Ti due to the opening 
of a gap. This was experimentally observed at the tran- 
sition between the IC high field phase and the uniform 
phase in the spin-Peierls compound CuGeC>3 |21|. Due 
to the IC nature of the ground state, the relaxation rate 
below T c should be dominated by the phasons, which are 
the standard Goldstone modes of IC phases, so that the 
decrease of 1/Ti is not necessarily thermally activated. 
This could be the origin of the apparent power law ob- 
served here. As shown in Fig. |J, a faster decrease is 
observed for H — 10.85 T, which corresponds approxi- 
mately to -f/cff = (H c i + H C 2)/2 where commensurability 
occurs. For H = 7.65 T, a field close to H&, the de- 
crease is noticeably slower. However, for this field value, 
1/Ti was measured on the line I [protons H(2) and H(4)], 
while it was measured on the line II for H = 8.0 T and 
10.85 T (see inset to Fig. 3). As discussed in Ref. [pL the 
corresponding sites have different form factors, and do 
not probe the same linear combination of the transverse 
and longitudinal spin-spin correlation functions. So, we 
cannot really attribute the weaker T-dependence of 1/Ti 
at 7.65 T to the proximity to H c \. 

In summary, for approximately the same T-range 
where 1/Ti decreases, we observe displacements of pro- 
tons located in the hydrogen bonding along the exchange 
path J||. These three features, namely i) absence of 1/Ti 
divergence at T c , ii) decrease of 1/Ti above T c and iii) ev- 
idence for proton displacements in the same T range, rule 
out any model involving solely a magnetic coupling be- 
tween the ladders Jt],|io| , and strongly support a " gener- 
alized spin-Peierls" scenario. We also found that the field 
dependence of T c is consistent with the predictions of a 
Bose condensation type of transition (T c oc (H—H c i) 2 / 3 ). 

It must be stressed that NMR spectra only tell us 
about the time-averaged displacements of these protons. 
Would the displacements of protons H(2) and H(4) be 
purely static, only the value of Jii would change (and 
thus that of H c i cx J± — J\\). The dynamics, evidenced 
by the divergence of 1/Ti in the ID regime, would not 
be affected. To alter the magnetic excitation spectrum, 
a dynamical modulation of the position of these protons 
must be present. In other words, they have to partici- 
pate to some phonon mode coupled to the magnetic ex- 
citations and leading to a dynamic modulation of Ju . 
This magneto-elastic coupling appears prior to the " spin- 
Peierls transition" at T C (H) and it readily explains the 
change in the T-dependence of 1/Ti above T c . The freez- 
ing of this collective mode would finally lead to a static 
IC modulation of the position of the protons along the 
legs of the ladder. 

We note that preliminary data show that anomalous 



shifts, related to proton movements, are still observed 
very close to H c \, where T c — > 0. This can be explained 
by the fact that, already at H = H c \, the system has in- 
terest in suppressing the quantum magnetic fluctuations 
p2[ through the magneto-elastic coupling. However, an 
ordering cannot be achieved at finite temperature since 
there is no finite magnetization along the individual lad- 
ders. 

We thank Th. Giamarchi for discussions and P. Van 
der Linden for technical help. 
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